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ABSTRACT 

The QCD corrections to the top-quark pair production via both polarized 
and unpolarized gluon fusion in pp collisions are calculated in the Minimal 
Supersymmetric Model(MSSM). We find the MSSM QCD corrections can 
reach 4% and may be observable in future precise experiments. Furthermore, 
we studied the CP violation in the MSSM, our results show that the CP vi- 
olating parameter is sensitive to the masses of SUSY particles (It becomes 
zero, when the cm. energy is less than twice the masses of both gluino and 
stop quarks.) and may reach 10 -3 . 
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I. Introduction 

The minimal supersymmetric model(MSSM) jl| is one of the most interesting ex- 
tensions of the Standard Model (SM). Therefore testing the MSSM has attracted much 
interest. As is well known, the MSSM predicts supersymmetric(SUSY) partners to all 
particles expected by the SM, and searching for their existence is very important. 

Since the top-quark was already found experimentally by the CDF and DO Collab- 
orations at Fermilab we believe that more and more experimental events including 
top-quark will be collected in future experiments. That gives us a good chance to study 
the physics in top-quark pair production from pp or pp collisions with more precise exper- 
imental results. Because of the heavy mass of the top quark this process provides a test 
of the SM and possible signals of new physics at high energy. 

The dominant subprocesses of top-quark pair production in pp or pp colliders are 
quark-antiquark annihilation and gluon-gluon fusion. The lowest order of those two sub- 
processes has been studied in Ref. ||. There it was found that the former subprocess 
(qq annihilation) is more dominant in pp collisions when the cm. energy(y / s) is near the 
threshold value 2mt, whereas subprocess via gg fusion will be more and more important 
with increasing cm. energy, and can become the most dominant one when the cm. energy 
is much larger than 2mt- 

In Ref. Q, the QCD corrections to top-quark pair production in pp collisions have 
been studied in the frame of the SM. It may seem natural that the QCD corrections of 
those processes in the frame of the MSSM are important for distinguishing those two 
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models. Recently, the SUSY QCD corrections to top pair production via qq annihilation 
were given in Ref. The SUSY QCD corrections via unpolarized gluon-gluon fusion 
were presented by C.S.Li, et. al ||. 

It is obvious that the correction from the SUSY QCD is related to the masses of top- 
quark and SUSY particles. Assuming the SUSY breaking scale at about 1 TeV, the masses 
of SUSY particles would be smaller than 1 TeV. So we can hope that corrections from SUSY 
particles are significant, since the heavy mass of the top quark (mt = 175.6 ± 5.5 GeV 
(world average)) may be comparable to some of the light SUSY particle masses. Therefore 
the SUSY QCD correction would give us some significant information about the existence 
of SUSY particles indirectly. 

Recently, the spin structure of the nucleon has been intensely studied by polarized deep 
inelastic scattering experiments at CERN and SLAC. This knowledge allows us to find a 
clear signal beyond the SM, if we collect enough events in the process of top-quark pair 
production from polarized pp or pp collisions. In the SM QCD, there is no CP violation 
mechanism, whereas in the SUSY QCD, the situation may be different. If we introduce 
phase angle of quark SUSY partners, we can get CP violation in the MSSM QCD jjffl. 
Once we get enough statistics of top-quark pairs from pp or pp colliders at higher energy, 
it will be possible to test CP violation. On the other hand, the spin-dependent parton 
distributions can be obtained from their polarized structure function data in Ref. || [ 1C ] [ 1 1 ] . 
There one found that the shape of polarized gluon and quark distributions in the nucleon 
depends on its polarization. Therefore the CP violation effects through the process of 
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top-quark pair production via gg fusion may be observed in polarized pp or pp collisions. 

In this work we concentrate on the SUSY QCD corrections to the process pp —* gg — » 
tiX both in polarized and unpolarized colliding beams. In section 2, we give the tree level 
contribution to subprocess gg — * ti. In section 3 we give the analytical expressions of the 
SUSY QCD corrections to gg — ► ti. In section 4 the numerical results of the subprocess 
gg — * ti and the process pp — > gg — > ttX are presented. The conclusion is given in section 
5 and some details of the expressions are listed in the appendix. 

II. The Tree-Level Subprocess 

The graphical representation of the process g(Xi, k{)g{\ 2 , k 2 ) — ► t(pi)i{p2) is shown in 
Fig.l (a). The Mandelstam variables are defined as usual 

a =(pi+P2? = (h + k 2 ) 2 (2.1) 
i ={pi-h) 2 = {k 2 -p 2 ) 2 (2.2) 
u = (p x - k 2 ) 2 = (h - P2 ) 2 (2.3) 

so s + i + u = 2m 2 . The amplitude of tree- level diagrams with polarized gluons can be 
written as:Q (a, b are color indices of external gluons , i,j are colors of external top-quarks 
and T a = 4^ are the Gell-Mann matrices.) 

M W = gy^uh)e v >\\2,k2)u i (p 1 )r( i y Vj {p 2 ), (l = s,t,u) (2.4) 

with 

r W = _ h)g ^ + (2k 2 + h)^ - (2h + k 2 ) vlp \ (2.5) 
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r(t) = _EA lii{h _ h + mt)lv (2.6) 
r( „) = lE^ lv{h _ h + mt) ^ (2 .7) 
We chose a form in which only physical polarizations of gluons remained: 

e^^hVi^h) = %*(-<r + ^ K X K - tSf + iX ^ pv (2 - 8) 

where n = fci + &2, Ai^ = ±1. /.From that, we can get the cross section at the tree-level 
with both polarized and unpolarized gluons. 

III. SUSY QCD corrections (non-SM) to the subprocess gg — ► ti 

1. Relevant Lagrangian in the MSSM. 

The difference between the MSSM QCD and the SM QCD corrections stems from the 
interactions of SUSY particles. Thus we can divide SUSY QCD corrections into a standard 
and a non-standard part. The Lagrangian density of the non-SM part of the SUSY QCD 
interaction is written as: 

L =Li+L 2 + L 3 + L 4 (3.a.l) 

Where 

L x = -ig s A»Tf k {q{d^ k L ~ <?1^<?D + (L - R) (3.a.2) 
L 2 = -y/2g s Tf k {l a P L q k q j * + <?P R g a q\ - ~g a P R q k q R * - q> ' PlMr) (3.a.3) 

L 3 = ^g s fabJ a ^9 b A^ (3.a.4) 
U = y 2 s A«A%(q* L q L + q* R q R ) + y^ c A^ b (qfT^ + $^4) (3.a.5) 
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q stands for quark, q for corresponding squark, g for gluino, Pl and Pr for left, right 
helicity projections, respectively. The mixing between the left- and right-handed stop 
quarks tL and Ir can be very large due to the large mass of the top quark, and the lightest 
scalar top-quark mass eigenstate t\ can be much lighter than the top-quark and all the 
scalar partners of the light quarks. Therefore the left-right mixing for the SUSY partners 
of the top quark plays an important role. Here we only considered the SUSY QCD effect 
from stop-quark, because we assume that other scalar SUSY quarks are much heavier than 
the stop-quark and hence decoupled. Furthermore we introduce the phase angle <pA in the 
stop mixing matrix. Defining 9 as mixing angle of stop-quark, we have 

t L = e^(£i cos 6> + t 2 sin #) (3.a.6) 

t R = e~ (-tisin6» + t 2 cos0) (3.a.7) 

where we suppose < . 

2. Analytical results of the MSSM QCD corrections. 

The one- loop SUSY QCD correction diagrams are shown in Fig. 1(b). In the following 
we present only the amplitude expressions of s-channel and t-channel. The amplitude of u- 
channel can be obtained from the t-channel expression by the following variable exchanges: 
t <-> u, k\ <-> foj e^(fci) <-> e^(/c2) and T a <-> T b . The one-loop diagrams can be divided into 
three groups: the self-energy diagrams of gluon and top-quark shown in Fig.l(b.l); gtt and 
ggg vertex correction diagrams shown in Fig.l(b.2); box diagrams shown in Fig.l (b.3). 
The ultraviolet divergence is controlled by dimensional regularization (n = 4 — e). The 
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strong coupling-constants are renormalized by using the modified Minimal Subtraction 
(MS) scheme at charge-renormalization scale /i/j. This scheme violates SUSY explicitly 
and the qqg Yukawa coupling g s , which should be the same with the qqg gauge coupling 
g s in supersymmetry, takes a finite shift at one-loop order. Therefore we take this shift 
between g s and g s as shown in Eq.(3.b.l) into account in our calculation, in order to have 
the physical amplitudes independent of the renormalization scheme and we subtract the 
contribution of the false, non-supersymmetric degrees of freedom (also called e scalars) 







g s = g s 



4vr 3 2 ' 



(3.6.1) 



where Ca = 3 and Cp = 4/3 are the Casimir invariants of SU(3) gauge group. The 
heavy particles(top quarks, gluino, stop-quarks, etc.) are removed from the jiR evolu- 
tion of a s (/i|j), then they are decoupled smoothly when momenta are smaller than their 
masses [13]. We define masses of heavy particles as pole masses. 



The renormalized amplitude corresponding to all SUSY QCD one-loop corrections (as 
shown in Fig.l) can be split into the following components: 

5M = SM S + SM V + 5M box + 5M d . (3.6.2) 

where 5M S , 5M V , dM^ox and 5Md are the one- loop amplitudes corresponding to the self- 
energy, vertex, box correction diagrams and the decoupling part, respectively. The 5Md 
stems from the decoupling of the heavy flavors from the running strong coupling, and is 
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given explicately by (see also |12| 

5M d = M (^M)[i log(£f ) + i log(£* ) + +i log(^) + | log(§)] (3.6.3) 
3. Self-energy corrections to the amplitude. 

The amplitude of self-energy diagrams 5M S (Fig.l.(b.l)) can be decomposed into <5Mf 

(gluon self-energy) and 5Mj (top-quark self-energy), i.e. 

5 M s = 5M§ + 5M§ 

= 5Ml {s) + 8M°® + 5Mf u) +SMf' ) +5Mf u) . (3,&1) 

The amplitudes 5M 9 s (s \ 5M 9 s {t) and 5M 9 s (u) are for s-, t- and u-channel, respectively. They 

can be expressed as: 

5Ml {s) = iM (s) [n(^) + Ii(k 2 2 ) + 2II(s)], (3.C.2) 
SMf (t) = ±M (t) [n(& 2 ) + u(k 2 2 )], (3.c3) 
6Mf u) = ^ u) [U(kj) + life 2 )], (3.c4) 

where Mo is the tree- level amplitude defined in Eq (2.4). 

n(fc 2 ) = -^(T F (B + AB 1 +4B 21 )[k,m il ,m il ] + 

T F (B + 4B 1 + AB 21 )[k,m i2 ,m- t2 ]-4C A (B 1 + B 21 )[k,rng,m~ g } - \C A ). 

(3.C.5) 

where Cp = I, Tp = \ , Ca = 3 are invariants in the SU(3) color group, Bi and Bij 



are Passarino-Veltman two-point functions [|J]|15]. The definitions of Bq, B\ and B 2 \ are 



listed in Appendix A. The amplitude 5Ms^ is written as: 

6M!® = --^^e^y^u^, (3c6) 

(fa-fa + m) ^ki(k2-P2) (h-fo + m thvVj(p 2 ). 



Here we define 



Sfci(p) = C F {H L {,P L + H R {>P R - H s l P l - H s R P R )5 kl (3.c7) 



with 



"2 p 

#L = ■§ftx 1 xzB 1 \p,mg,m il ] + (m ?1 -> m t ~ 2 ,Xj -> j/j) + ±(5Z L + &Z£), (3.C.8) 

#R = ^^2^4-Bi[p, mj, m t -J + (m fi -> m^x* -> j/i) + t^Zr + (3.c.9) 

fff = fi^sm^m^J + K x - m h , Xi - w ) + + «4) + <5m t) 

(3.C.10) 

H R = J^^i x 4"T-g-Bob, m t -J + (m^ -> m t ~ 2 , Xj -»■ yj) + \m t {5Z R + <5z£) + <5m t . 

(3.C.11) 

Where we abbreviate (f> = 4>a, x\ = cosfle -1 ^, X2 = sinOe 1 ^, x 3 = cosOe 1 ^, x 4 = sm9e~ 1 ^, 
yi = sin^e - *^, 1/2 = —cosOe 1 ^, y 3 = sin#e^, y 4 = — cosde~ l< ^, and is mixing angle of 
stop-quarks, see Eq (3.a.6 ~ 7). 

The explicit expressions of the top-quark wave- function renormalization constants have 
the following forms: 

-2 

5Z L = -^7(xix 3 i2e[Si] - ^j(xix 4 - x 2 x 3 )Re[B ] 

+ m t 2 (xix 3 + x 2 x 4 )-Re[ J B 1 ] (3.C.12) 



m t m g (x 2 x 3 + xiX4)Re[B' ])\p,mg,m {l ]\jp =m 2, 



(3.C.13) 



8Z R = -£^(x2XiRe[Bi] + m\{x\x 3 + X2X4)Re[B' 1 ] 
- m t m g (x2X 3 + XiX 4 )i?e[ J Bo])[p,mg,m t ~ 1 ]| p2=m 2, 

"2 

5mt = xf^((xiX3 + X2X i )m t Re[Bi] (3 c 14) 

- (x 2 x 3 + xix 4 )mpi?e[ J Bo])[p,mg,m t ~ 1 ]| p 2 =m 2, 



We use the following abbreviations: B i i ^[p,m\,m2\ = 9Bi ^^P^ ll ' m2 ^ _ 
4.Vertex-corrections to the amplitude. 
The amplitudes for vertex diagrams can be expressed as: 

6M® = g s e^ a {k 1 )e^ b (k 2 )u i (p 1 )Amv j (p 2 ), (I = s,t,u), 



(3.ci.l) 



where 



A(') 



fa 



&{k u k 2 ) 



7p 



r[(h - k 2 ) p g^ v + {2k 2 + ki)^g vp 



-{2ki + k 2 ) v g t ip\ 



A' 



P,(ij) 



(Pl,P2y 



and 



AW = t^k {^mj [A« ( im) (Pl, h ~ Pi)} (h ~h + TTHhu 
+T? m ^(h ~i> 2 + m t ) [K,(mj)( k 2 - P2,P2)] } ■ 



(3.d.2) 



(3.d.3) 



The functions A^ and A^ ^ are listed in Appendix B. 
5. Box-corrections to the amplitude. 

The box diagram corrections in the t-channel (Fig.l(b.3)) are given as follows: 
SM®, = 2g^ a (k l )e^\k 2 )u i {p l )({T c T a T b T% 3 F { ^ ) 



ifbcd{T c T a T^)ijF^ u ^ facmfbmd{T c T^')ijFj] lU ^ 

-[T c (T a T b + T 6 T a )T%-F^ 4) K(P2), 



(3.e.l) 



Where f a i, c is defined as [T a ,T b ] = if a bcT c - The form factors F$(i = 1 — 4) correspond to 
the kernel of the four Feynman diagrams in Fig.l(b.3) respectively and are given explicitly 
in Appendix C. 



6. Total cross section. 
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Collecting all terms in Eq (3.b.2), we can get the total cross section: 

<7(Ai,A 2 ) =a (A 1 ,A 2 )(l + ^(A 1 ,A 2 )) 

= £ dtZ spins [\M \* + 2Re{Ml5M)] [ ^ ' 

where = {ml — \s) ± \sj3t, fit = \A — 4mf /s, and the spin sum is performed only over 
the final top-quark pair when we considered polarized gluons. 

IV. Numerical results 

We denote gq for the Born cross section and a for the cross section including one- 
loop SUSY QCD corrections of subprocess gg — ► ti, and define its relative correction as 

6 = ■ For polarized gluon fusions, <r++, a and <r+_ are the cross sections with 

positive, negative and mixed polarization of the gluons, respectively. In order to inspect 
the CP violating effects we introduce the CP-violation parameter for the subprocess defined 
by icp = a^+all • The possible SUSY QCD effects in gg — ► ti should be observed in pp 
colliders. By analogy we can define also the relative correction and CP violating parameter 
for the process pp — > gg — > ti as S = a ~°° and £cp = a+++a ' respectively. The SUSY 
QCD contribution to the process p(Pi, x)p(P 2 , y) — ► gg — > iOf (x,y are polarizations of 
protons) can be obtained by convoluting the subprocess with gluon distribution functions. 

a(s) = J dxidx2G(xi,Q)G(x2,Q)a(s,a s (iJ,)) (4.1) 

with k\ = x±Pi, fc 2 = x 2 -P 2 and r = xix 2 = s/s. G(xi,Q)(i = 1,2) are gluon distribution 
functions of protons. We take Q = [Xr = 2mt- 
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In order to get results of top quark pair production from polarized pp collisions, we need 
to consider the polarized gluon distributions in protons. The cross sections of polarized 
PP ~* 99 ~^ ttX can be written as 

a(x,y) = Y,x 1 ,\ 2 =± J dx 1 dx 2 G xXl (xi,Q)G yX2 (x2, Q)&\ lt \ 2 (s,a s (p>)) (4.2) 

where x and y are the polarizations of incoming protons and Ai, A2 are the polarizations 
of gluons inside protons. G xXl (x,Q), G yX2 (x,Q) = G ± (x,Q) for equal (+) and opposite 
(-) polarization, G + (x,Q) and G~(x,Q) are polarized gluon distribution functions in the 
proton. 

We used unpolarized proton structure functions of Gliick et al. || in our numerical 
calculations. For the polarized proton structure functions, we use the evolution equations 



of Gliick et al. [10| with input parameters from the paper of Stratmann et al. JII] (Next- 
To-Leading-Order). Since the structure functions are one of the least certain input in our 
calculation, we checked the result against other set, i.e. the polarized structure functions 
G ± (x, Q) of Brodsky et. al (8| (Using Leading-Order only). This tests the stability of our 
results against the particular form of the input structure functions. The two different sets 
of input are compared in Fig. 2, which gives the relative SUSY QCD correction {5) and 
^CP versus cm. energy for the process pp — > gg — ► tiX. Though the SUSY QCD 
corrections from the two sets of structure functions are not too different, for 6, there is 
some noticable change for £cp- Because ^cp depends strongly on the cm. energy of 
the subprocess gg — * tt (shown in Fig. 3(b)), a small modification of structure functions 
may lead to a large change of £,cp- Thus we can infer that the NLO-QCD calculation 
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is required and the precise numerical prediction does depend on the reliability of the 
structure functions. 

The SUSY QCD relative corrections are about 2% ~ 4% and decrease with increasing 
cm. energy from Fig.2. These correction effects are within reach of future precision 
experiments and provide a possible discrimination of the SM and the MSSM effects. From 
Fig.2(c) we can see that the CP violation parameter ^qp can be 10~ 3 . Therefore, CP 
violation in this process stemming from the SUSY QCD can in principle be tested in 
future precision experiments. That would help us to learn more about the sources of CP 
violation. 

In order to explore the effects of the SUSY QCD correction for future arrangements of 
optimal experimental conditions, we also investigate the subprocess gg — ► ti. 

The relative SUSY QCD correction and CP violating parameter versus cm. energy 
(V$) for different polarization gluons are plotted in Fig.3 (a ~ c) with m g = 200 GeV, 
m h = 250 GeV, m h = 450 GeV, and 6 = = 45°. In Fig.3(a) 5 ++ and <L_ 
are drawn in solid line and dashed line, respectively, ^cp as a function of cm. energy 

is depicted in Fig.3(b) and 5^ as function of VI is plotted in Fig.3(c). Each curve in 

Fig.3(a) has an obvious peak near the position of the threshold of top pair production. 
That large enhancement is the combined effect of the threshold, when a/I is just larger 
than 2m t = 350 GeV, and the resonance when \f§ ~ 2rrig = 400 GeV. The small spikes 
around the position of \/i = 900 GeV, there V§ ~ 2m t ~ 2 = 900 GeV, shows also the 
resonance effect. Although Fig.3(a) shows that (5++ and 5 approach equal values when 
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the cm. energy is far beyond its threshold value 2m t , the quantitative difference between 

<5_| and <5 ++ still exists in the whole energy range plotted in these figures. Fig.3(b) 

shows also that £cp will be zero, if the cm. energy is below the threshold of SUSY 
particles in the loop( i.e. a/I < 2m g = 400 GeV in Fig.3(b)). This is reasonable 
because only beyond this point can we have absorptive terms which give contributions 
t° S.CP- icp has obvious resonance effect in the regions around \/I ~ 2m g = 400GeV 
and a/I ~ 2m t ~., (i = 1,2) = 500 GeV, 900 GeV. We also find that the two stop quarks 
give opposite contributions to £cp and when their masses are degenerate £cp will vanish. 
When the cm. energy a/I is larger than 1 TeV, ^cp will be near zero, because the 
contributions from the two stop quarks will cancel each other. Therefore a quantitative 
strong change of £cp as function of cm. energy can be an indication for the signals of 
stop quarks and gluino. 

<t(±, ±) and £cp as functions of rrig are shown in Fig.4 (a) and Fig.4 (b), respectively. 
In Fig.4 we take a/I = 500GeV, m h = lOOGeF, m t ~ 2 = 450G e y, 9 = <j> = 45°. We can 
see from Fig.4(b) that £cp changes its sign when m g is near m t = 175 GeV . The curves 
in Fig.4(a)(b) show again the resonance effect when a/I ~ 2m g = 500 GeV, note that for 
each line there is a steep change of the value of <r(±,±) or £cp around the position of 
m g = 250GeV. 

Dependences of relative correction 6±± and £cp for the subprocess gg — ► ti on 
are plotted in Fig.5 (a) and Fig.5 (b). 6±± and £cp as functions of are shown in 
Fig.6 (a) and Fig.6 (b), respectively. In all figures of Fig.5 and Fig.6, we take the common 
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parameter set with y/§ = 500 GeV, rrig = 200 GeV and 9 = <fi = 45°. In Fig.5, we set 
mi 2 = 450 GeV , whereas = 100 GeV in Fig. 6. We find that £cp m fact increases 
with mass splitting of stop-quarks (i.e. m^ 2 — m^) and when = m ( - 2 , ^cp is equal to 
zero. The resonance effect of stop quarks, when y/§ ~ 2 m,f (i = 1,2), is superimposed 
on the curves in Fig.5 (a)(b) and Fig. 6(a) (b) around the positions of = 250 GeU 
in Fig.5(a)(b) and = 250 Gey in Fig.6(a)(b), respectively. Around those points the 
relatively sharp changes of the values of £cp and the relative corrections are shown in 
these figures. 

Finally, the dependence of 6±± and £cp on the phase (f> is shown in Fig.7 (a) and (b). 
In Fig.7, we take yfl = 500 GeV, m g = 200Gey, 9 = 45° and m h = 150 GeV. We find 
that £cp is directly proportional to sin (20) and reaches its maximal value when 4> = \. 

IV. Conclusion 

In this work we have studied the one-loop super symmetric QCD corrections to the sub- 
process gg — ► ti and process pp — ► gg — > ttX. The calculations show that the SUSY QCD 
effects are significant. The absolute values of the corrections are about 2% ~ 4%, so they 
may be observable in future precision experiments. Furthermore, we find £cp depends 
strongly on masses of SUSY particles and can reach 10~ 3 when we take plausible SUSY 
parameters. 

The results show that there is an obvious difference between the corrections for the 
protons polarized with parallel spin and that with anti-parallel spin. Hence there is a 
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possibility to study spin-dependence in the frame of the MSSM QCD. 

We also presented and discussed the results of the subprocess gg — > tt. We find that 
when the cm. energy passes through the value 2m g or 2m^ (i = 1,2), the value of the 
CP violating parameter £cp changes strongly. If cm. energy is less than both 2m g and 
2rrif. (i = 1, 2), £cp will be zero. If in future experiments a sharp change in ^qp is found 
with \/i running from low cm. energy to high cm. energy, it would be interpreted 
as a signal of SUSY particles. Furthermore, because the CP violating parameter £cp is 
sensitive on the mass of gluino(as shown in Fig.4 (b)) and the mass splitting of stop-quarks 
mi 2 — mi x (as shown in Fig.5 and Fig.6), we can also get information of SUSY particles 
from precise measurements of £c*p. 

The authors would like to thank Prof. A. Bartl for useful discussions and comments. 
One of the authors, Yu Zeng-Hui, would like to thank Prof. H. Stremnitzer for his help. 
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Appendix 

A. Loop integrals: 

We adopt the definitions of two-, three- and four-point one-loop Passarino-Veltman 
integral functions of reference! 14] [15]. 



l.The two-point integrals are: 



{B ;B^B^ u }(p,mi,m 2 ) = / d n qy^ ' /" f st, (A.a.l) 

itt z J [q z — mf\ [(q + p) l — m^J 

The function £?„ is proportional to p^: 

Bfj.(p,mi,m2) =p fi B 1 (p,m 1 ,m 2 ) (A.a.2) 

Similarly we define: 

Bp* = PnPvB 2 i + g^B 2 2 (A.a.S) 

We denote B = B - A, B 1 =B 1 + ±A and B 21 = B 21 - ±A. with A = f - 7 + log(4vr), 
e = 4 — n. is the scale parameter. 
2. The three-point integrals are: 

{C ;C fM ;C f _ lu ;C^ up }{p,k,mi,m 2 ,m 3 ) = 

7 [q 2 — mf][(q + p) 2 — m|][(g +p + A;) 2 — m|] ' 
We define form- factors as follows: 
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Cf,u = PpPuC 21 + k lx k v C 2 2 + {PpK + k^p^C^ + g^C 2 4 
Cp Vp = PpPvPpCzx + kpk v k p C 32 + (k^p u p p + p p k u p p + p p p v k p )C 33 + 
{kpKpp + ppk v k p + kfj,p v k p )C34 + {Pp9u P + PvQpp + P P 9pu)C35+ 
{kp9v P + Kg pp + kpg^Cm (A.a.5) 
3. The four-point integrals are: 

{-Do; Dp; D pi/ ; D^p; Dn„ pa }(p, k, I, m 1 ,m 2 , m 3 , m 4 ) = 
(27r^) 4 ~ w /■ {1; g M ; q p q y ; qpq u q p ; g^gpfe} , . 

ivr 2 7 (? [ g 2_ m 2 ][(g+;p) 2_ m 2 ][((7+;p+A , )2 _ m 2 ][(g+p + A , + /) 2_ m 2 ] ' I 

Again we define form- factors of D functions: 

Dp = PpDu + k p D 12 + IpDiz 

Dpv = PpPuD 2 i + k p k u D 22 + l p l v D 23 + {pk}p, v D 2i + {p/}^D 2 5 + {kl}pvD 26 + g pu D 27 
Dpup = PpPnuP p D 31 + k p k nu k p D 32 + lpl nu lpD 33 + {fcpp} Mi , p -D 3 4+ 
{ipp}pu P D 35 + {pkk} pup D 36 + {pll}p U pD 37 + {/fcfe} /ll/p D 38 + 

{fe^}^P^39 + {p^j^pAilO + {^1^^311 + {%}^z/p-D 3 l2 + {Zp}>z/p-C>313 (Aa.7) 
where 

{pk}pu = Ppk u + kpp v 
{pkV\ pup = Ppk v lp + IpPvkp + k p l v pp 

{pg}p,vp = Ppgvp + JVSV + PpSV (Aa.8) 
17 



The numerical calculation of the vector and tensor loop integral functions can be traced 



back to the four scalar loop integrals Aq, Bq, Cq and Dq in Ref.[14|[15| and the references 
therein. 

B. Vertex corrections: 

The 3-gluon- vertex can be written as: (a,b,c are the color indices of the external gluons) 



AiS(fcl, fea) = j£z [Tr{T b T c T a ) [k%{k u fc 2 )] + if<™n f anl f blm [ A g>p(fci, k 2 )] } 



(A.b.l) 



the vertex functions A^p, A^] p are expressed as follows: 



A$p(fci, k 2 ) = fig^kiv + f^g^hp + fi gu P k 2lx + fi a) g^k 2p 



f( a ), 



r(a). 



(A6.2) 



where a = 1,2, and the /i 1 ^, /i 2 ^ are given in terms of the Passarino-Veltman functions with 



internal stop lines C-j (= Cjj[— — fe, m^, rrifj) and internal gluino lines C, 
Cjj[— fci, — /s2 5 m g-> m gX)- For simplicity, we abbreviate the definite part of C integral 
functions (using the definitions of @ ||) as follows: C$ = C$ - ±A, C^f = C$ + ±A, 



(2)/ 



^?=C^ + ^A. (a = 1,2) 



(i) 



-8(^24^ — 8C' 



(i) 

35 ) 



r(l) 



-8a 



(i) 



36 > 



(A.b.3) 



J4 - _4U 24 - °°36 ) 
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#=4C« + 12C« + 8Cg\ 



/ 6 W=4Cg)+8Cg) + 4Cg) + 8Cg), 



and 



A {2) = -8m?cf - 4m|cg> - 16c£> + l2eC$ - 8C^ + &cg> 
+8h ■ k 2 C$ + I6&1 • + 8/ci • k 2 C$, 

fW = -impif - 8C^ + 6eC^ + 8C$h ■ k 2 + 8C^h ■ k 2 , 

ft = 4m|^ 2) - 4m|cg> + 8Cff - 6cC® - 8Cg> + &(?£> ft 
+8A; 1 -fe 2 cg ) + 8fei-A:2Cg ) , 

/f = -4m§cf - 4m§cg> - 8C£> + &c£> - 8Cg> + fccg'4- 
+8^ • A; 2 C{ 2) + 8fci • k 2 C { £ + 8fci • k 2 C$ + 8fci • feC^, 

ft = "8Ci 2) - 24Cg } - 16C$ 

Q/^(2) Q/^P) 1 

J 6 — — »0 12 — 10O 22 — 80 23 - l0O 34 

Similarly, the gtt vertex functions are composed of left-handed and right-handed con- 
tributions plus a counterterm: (We define a as color index of external gluon and i,j as 
colors of external top-quarks) 

A« (ij) (pi,P2) = -§£T%{(2C F - C A )(A£ L \ Pl , P2 )P L + A£ R) ( Pl ,p 2 )P R ) 

+C A (A^ L \ Pl , P2 )P L + A^ R \ Pl , P2 )P R )} (A6.5) 
+(m^ -> m t - 2 , a;, -> yi) + A^f T) 

The expressions of A{ t n) ,n= 1L,1R,2L,2R are given as following: 
We define 

Co 3) , cff = Co , Cj \- P i , -p 2 , m fl , m-g , m ?1 ] 
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Cf ,cf = Co,C l3 [- Pl ,-p 2 ,m s , mii ,m § } 
Then we can get as follows: (i=l,2,..., 10) 

/if L) = -2x 2 x A c!$ 



4 1L) =x 2 x 3 mg(C^ + 2Cg } ) 



and 



4 1L) =x 2 X3m 5 (ci 3) + 2C^ } ) 
4 1L) = ^(cf } + 3Cf? + 2Cg>) (A.6.7) 
4 1L ) = x 2 x 4 (cf } + c[f + 2(7$ + 2Cif) 
h£ L) =x 2 x 4 (cif + 2Cg>) 
4 1L ) =x 2 x 4 (C'g ) + 2Cg ) ) 

4 1L) = 4 1L) = 4o L) = o 



h[ 2L) = x 2 x,(-mjC^-2C^ + eC^) 

+x 2 x 4P j(c£ ) + C$) + 2x 2 x m p 2 (C^ + 4?) 
+x 2 x,pl{C i S + 



hf L) = 2x 2 x 3 mjc{t ] 
hf L) = -2x 2 x 4 (c{f + C^) (A.b.8) 
4 2i) = -2x 2 x 4 (cfJ + cg>) 
4 2L) = -2x 2 x 4 (ci? + cg ) ) 
20 



4 2L) = -2x 2 x 4 (c{J + d§) 



= = x^m\cf 

/4 2 L) = x 2X Mt ] - C[f) 

h^ R ^ and h\ 2R ^ can be obtained by exchanging x\ <-> a?2 and X3 <-> X4 in h\ 1L ^ and /ij 2i ' ) . 
(i = l,2,..,10) 

The counter terms are given by: 

kf T) = -CpSfT^ [{5Z L + 5Z{)P L + (SZ R + 8Z R )P R ] (A6.10) 

The wave function renormalization constants can be obtained from Eq.(3.c.l2) and Eq.(3.c.l3). 
C. Box corrections: 

Finally, we list the four form factors Fjf v as given in Eq.(3.e.l) in terms of Passarino- 
Veltmann functions. First, we define F l k iL and F^ R by: 

+ P ^A tlR) +P2,iA tlR) 

+7,hPiAi R) + 7,hP2uF^ R) (A.c.l) 
+luhPlv F i3 R:} + liyhP2Ai lR) 

+PltA,P2uF^ R) + Pl^PluF^ 
+P2nPlyFfj R) + hPlvPlisFi 8 lR) 
+hP2^P2uF^ R) + P2 l xP2 V F^ R) ] 

+{Pr - Pl, F^ tiR) - if L) ), (k = 1 ~ 20, i = 1 ~ 4). 
In the following we only give the expressions of Fl iR (k = 1, 2, 20 and i = 1 ~ 4). The 
expressions of Fj? L can be obtained from F^ R by exchanging x\ <-> X2 and X3 <-> X4. 
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Furthermore, the form factors in the u-channel are given by 



F^ ) (k 1 ,k2,Pi,P2) = Fff(k2,k 1 , Pl ,p 2 ) (A.c.2) 



The expressions of F^ tlR ^ (k = 1 ~ 20) are given as below: 

F (tiR) = p (tiR) 



= —2xiX4mgD^ + 2x 2 X4m t D^ 1 (A.c.3) 
+ 2(xix 3 - x 2 X4)m t D^ 3 , 

Ft R) = ^(Dffi-Dg)), 
Ft R) = -4x lX3 (D^ + D^ 2 ), 
Ft R) = 4x 1 x 3 (^ + <-^ 3 ), 
Ft R) = -4x^3^1 

F (u R) = p(tiR) = 2xiX3{D m _ D w )t 

F^ = 4x lX3 (D$ - D<J - D% - D« - D% 

+DU + 2DU-D® + DV-Dj8 + DU), 

F£ R) = 4x iX3 (D^-DS + D^-DS ), 
p(tlfi) _ „(tlR) _ p (tlR) _ UtlR) _ n 

F^ = ^xsMDS - + 2>£> + ^ - i$ + Dg) ) 
-4x 1 x 4 m § ( J DS 1) + + - £>£ + D« - £>$) 
+4x 2 x 4 m t ( J Df 1 ) - £>$ + £>£ } + + £>$ - 2D^ 
-^ + ^ + ^-2^ ), 

i#*> = 4x 1 x 3 m t (-^ + Off - D® + D&> - D« + Z>£> ) 
+4x 1 x 4 m 5 ( J Dii ) - £>£> + - _ D W + ^J) 
-4x 2 x 4 m t ( J Dy ) - D$ - + £>« + 

^34 Ziv 35 + -^37 "^39 + Z1J 310)' 
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F$ R) = 4x 1 x 3 m t (Z)g ) - I?39 ) - Ax lX4 mg(D^ - D<J) 
+4* 2 x 4 m t (Ag - £>$ + - £>$,), 

^ = 4x 1 X 3 (-^ + ^ ) -^ 1 5 ) +^ 1 6 ) +4 1 4 ) 

-^£ ) -^ ) +^ ) +^S ) -^ ) ), 

F^ = 4x 1 x 3 (-^ 1 3 ) + ^ 1 6 ) -^ 1 9 ) + ^ 1 8 ) ), 

F 2 ( * 1R) = 4x 1 x 3 m t (-^ 1 3 ) - Ajg ) + 4x 1 x 4 m 5 ( J Di 1 3 ) + £>$) 
+4x 2 x 4 m t ( J D« - 4 ] + DS - 

where we denote £)f ) , D^f , = A, Aj, Ajfc [~Pi ,h,k2,m g , m t ~ , m t ~ , m t ~ ] . 
The expressions of F^ 2R) (jfe = 1 ~ 20) are as follows: 

= 2xix 3 m t (A?7 + Af 13 ) + 2xix 4 m 5 £)g ) - 2x 2 X4m t (A?7 + ^Si). ( Ac - 4 ) 
p(t2R) _ 2x 1 X3m t Z?3^3 + 2xix 4 mgZ?27^ — 2x2X 4 mj A31I' 
Ff fi) = 4x 1 x 3 (-A$-AS + A? 1 ) 2 ), 

Ff R) = 4x 1 x 3 (Ag 2 -A? 13 ), 
Ff R) = 8*1x3(0$ + Agi) + 2x 1 x 3 mJ(lf ) + A??) 

+2x 1 x 3 m?(-A! 2) " ^ " 24? " 0$ " ^ " ^? " 4?) 

+2(x 2 x 3 + xix 4 )mjra t (Dj 2) + £>g } ) - 2x 2 x 4 m 2 (r4 2 i ) - D<§ + Aj? - A25 ) 

+4xix 3 fci • Pi (A?? + 2Al? + A^?) 

+4x 1 x 3 fc 1 • p 2 (A! 2) + Ai? + D$H + £><?„) - Ax lX3Pl ■ P2 (D$ + 2D$> + 

Ff = 2x 1 x 3 m|0g ) 

-2xix 3 m t 2 (Ai 2) + A2? + Dig + Ajs) + 2(x 2 x 3 + x lX4 )mgm t D^ 
+2x 2 x 4 m t 2 (A^ 2) - Ai?) + 4ziX3fci -Pi(£>$ + A5 2) o) 
+4x1x3^! • p 2 (A2 3 + A39) " 4xix 3 pi • p 2 (A2 2) + AP7) 



+8xix 3 Ay[ 3 , 



Ff R) = 2x 1 x 3 (0g ) + 0fi ) 2 ), 
Ff = 2x 1 x 3 A5? 2 , 
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,{t2R) 
9 



2x 1 x 3 m t (-D ( $ + - D$ + Dg } ) + 2x 1 x 4 m § (Dg ) + £>g>) 
-2x 2 x 4 m t (D^ - D { $ + D { $ - D$), 

2x\x 3 m t {D^ - Dg)) + 2xiX4m~ g Df^ - 2x 2 x±m t {Df^ - Dg)), 

4x lX3 m t (D$ - D<§ - Dg> + Dg) - Dg> + flg - £>g J + flg) 
+4x 1 x 4 m § (Dg ) - £>g } + Dg } - £>g } ) 
-4x 2 x 4 m t (Dg) - D$ + 2L>g ) - 2£>g> + £>g> - Dg)), 

= 4x 1 x 3 m t (^g) - + Dg) - 2Dg> + D® - Dg> + 4? ) 
+4x 1 x 4 mg(-Dg ) - 2£>g> + £)g } - Dg) + L»g } ) 
+4x 2 x 4 m t (Dg ) - + 2D$> - 2L>g ) + Dg) - £>g>), 

Fg 2 R) = 4x 1 x 3 m t (-Dg ) + - £>$ + £>g>) + 4x 1 x 4 m 5 (-Z)g ) - £>g> + £>g>) 
+4x 2 x 4 m t ( J Dg ) - Dg) + ^ - £>g> ), 

F^ = 4x 1 x 3 (Dg)-Dg)-Dg) + Dg)), 

F^ = 4x 1 x 3 (-Dg) + Dg)+Dg)-Dg)), 

Fw R) = 4x 1 x 3 m t (-L»g ) + Dg) - £>g> + £>g>) + 4x 1 x 4 m § (-Dg ) + £>g>) 
+4x 2 x 4 m t (D^ - Dg) + Dg) - Dg) ), 

where Df ] , D^ 2) , D^f k = A , Aj , Ajfe [-f>i , &i , ~P2 , , m t ~ , m t ~ , nig] . 



? (t2R) 



13 



„(t2R) 
^14 

P (*2il) _ 
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The expressions for Fjf 3R ^ (k = 1 ~ 20) are written as: 



if 3K) = 2x lX3 m t (D^ + 2D®) + x^mtmf (Z^ 3) + D®" 



- Xl x 3 mUD^ + 2D}? - D\y + + + 2£) w _ 

+2x1X4/71,31)27 + xiX4m|Z?Q 

-xix 4 m t 2 mg(I)J 3) + - 2Z)g ) + D§> + 2D$ - 2D®) 

+2x 2 x A m t {D® + 2£)g ) 1 - 2£>g ) 3 ) + x 2 x 4 m t m|(Z)i 3 1 ) - 

-x 2 x 4 m t 3 ( J Df 1 ) - £>g> + 2£>£ ) + 2D® - AD® + - 2D® - 3D® + AD®) 
+2x 1 x 3 m t A; 1 • p^Dg - D® - D® + £>g> + D® - D® - D® + D®) 
+2x 1 x 4 m g k 1 ■ Pl (D® + D® - 2D® + D® + D® - D® - D® ) 
+2x 2 x 4 m t k x ■ P i(D® - D® + D® + 2D® 

+n® - m® - n® - n® + n® - n® + 2n® + n® - 2n® ) 

+-U 2 4 6JJ 25 1J 2Q u 33 + -^34 ^35 + Zi ^37 + ^39 Z1J 310) 

+2x l x 3 m t k l ■ P2 (-D® - D® + D® - D®) 
+2x 1 x 4 m~ g k 1 ■ p 2 (~D® + D® - D®) 
+2x 2 x A m t k 1 ■ p 2 (D® - D® - D® + D® + D® - D®) 
+2x x x 3 m tVl • p 2 {D® + D® - D® + D®) 
+2x 1 x A m gPl ■ p 2 (D® - D® + D®) 
+2x 2 x 4 m m ■ P2 (-D® + D® + D® + D® - D®), 



(A.c.5) 



p (t3R) = _ 2 XlXim ~D® - 2 Xl x 3 m t D® 3 - 2x 2 x 4 m t (D® + D® ± - D® 3 ) 
r (t3R) ...... 

^3 



— — z,xix 4 //tgj^27 ~~ ^^1^3"H^3i3 — ^x 2 x 4 iibt\i^ 21 ~r ^311 313/' 

= A Xl x 3 {D® + 2D® + D® - 3D®) + 2 Xl x 3 mj(-D® + D® - D®) 
+2x 1 x 3 m>{D® - D® + 2D® - D® + 2D® + 3D® - AD®) 
—2x 2 x 3 m t m g D® — 2xix 4 m t mgD® — 2x 2 x 4 mf(D® + D®) 
+Ax 1 x 3 k 1 ■ Pl (D® + D® - D® - D® 
-D® + D®-D® + 2D® + D®-2D®) 
+A Xl x 3 h ■ P2 (-D® + D® - D® + D® + D® - D®) 
+A Xl x 3Pl • p 2 (D® + D® - 2D®), 

Fp R) = A Xl x 3 D® 2 - 2 Xl x 3 mp® 

+2x l x 3 m1(D® + D® + 2D® - 2D®) 
—2x 2 x 3 m t m g D® — 2x\x 4 mtm g D® 
-2x 2 x A m 2 t (D® + D®>) + 4x^3^ • pi{-D® + D®) 
+Ax 1 x 3 k 1 ■ P2 (-D® + D®) + A Xl x 3Pl ■ p 2 (D® - D®), 



25 



if 3R) = 4x 1 x 3 (-D®+D®) + 2x 1 x 3 mjD® 

-2xix 3 m 2 (D® + D® ) + 2x 2 x 3 m t mg(D^ + D® - D®) 
+2xix A mtmg{Df ] + D® - D®) 
+2x 2 x 4 m 2 (D® + 2D® - D® + D® - D®) 
+4x 1 x 3 A; 1 -p 2 (D®-D®), 

if 3R) = -Sx^D® - 2x 1 x 3 mf,D® 

+2x!X 3 m 2 (D® + 2D® + D® - 2D®) - 2x 2 x 3 m t m § D® 
—2x\x A m t mgD^ — 2x 2 x 4 m 2 (D® + D®) 
+4x 1 x 3 k 1 • Pl (D® - D® - D® + D® + D® - D®) 
+Ax 1 x 3 k 1 ■ p 2 (-D® + D®) + 4x!X 3Pl • p 2 (D® - D®), 

b 1 - - 4xix 3 (^ 312 - V 313 ) + xix 3 m g (D Q ~L> l2 +U l3 ) 
+x 1 x 3 m t 2 (-D® + D® - flg) - D® 

+2D® - 2D® - 2D® + D® - D® + 2D® + 2D® - 2D®) 
+(x 2 x 3 + xix 4 )m t mgD® + x 2 x 4 m 2 (D® + D®) 

+2x l x 3 k l ■ Pl (-D® - D® + 2D® + D® - D® - D® + D® - 2D® + 2D®) 

+2x lX3 h • p 2 (D® + D® - 2D®) 

-\-2x\x 3 pi ■ p 2 (D® - D® - D® + D® + D® - D®), 

Ff R ) = 2x 1 x 3 (D® + D®-D®), 
p(*3«) _ 4 ™m® + D® - D® - D® + D® - D® + D® - D® 1 

^9 — 4XiX 3 (±/ 2 2 + -^23 -^25 ^26 + -^36 ^38 + -^39 -^310 A 

F (*3A) _ 4a , 1 ^m(3)_ jD (3)_ £) (3)_ D (3) ^(3) ^(3)^ 
r 10 ~~ ^1-^3^25 -^26 -^37 -^38 ^^39 ^-^310/' 

fff^ = 2xix 3 m t (-D® + D®) + 2x 1 x 4 m § (-D® + D®) 

+2x 2 x 4 m i (-D® + D® - D® + D® + D® - D®), 

F^ R) = 2xix 3 m t (D® +D®) + 2xix 4 mgD® 
+2x 2 x 4 m t (D® - D® + D® - D®), 

Fg 3ii) = 2x!X 3 m t (D® - D® + D® - D®) + 2xix 4 m 5 (Dg ) - D®) 
+2x 2 x 4 m t (D® - D® + D® - D® - D® + D®), 

F (t3R) = _ 2xiX3 m t {Df 3 ) + D®)-2x 1 x 4 m § D® 
+2x 2 x 4 m t (-D® + D®), 
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Ft? 1 * = -4 X1X3M D<$ - Dg> + D$ + J>g> - Dg) + Dg); 

+4x 1 x 4 m § (- J Df 1 ) - £>g> + £>$ - £>g> + £>g>) 

+4x 2 x 4 m t (- J Dg ) + L>g } - £>g } - £>g } 

-D (3) + 2i} (3) 4 D (3) - D {3) - D (3) + 2D (3) ) 
U 2 4 "I" -^25 + ^26 -^34 -^39 + ^310 J' 



^16 = -^m t {D^-D^ + D^-Dfi-Dfi + D^-D^ + D^) 
+Ax lXA m~ g {-D { $ 4 £>g> + D { S - £>g> - + L>g } ) 



+4x 2 x 4 m t (I)g ) - I?2? " ^2? + ^2? 
, n(3) _ r,(3) _ 2£) (3) , r>(3) _ ^(3) 2£) (3) , 



p(*3i?) , _ _ m (3) n (3) n (3) n (3)v 

+4x 1 x 4 m,( J Dg ) - Z)g } + £>$) + 4x 2 x 4 m t (-Z)g ) + £>g> - £>g> + £>g|,), 

^18 — 4X1X3 (,-^22 "^24 + -^25 -^26 -^34 + -^35 + -^36 
_ D (3) _ r)(3) D (3)x 
-^37 -^38 + -^39 /) 

r 19 ~ 4X1X3^23 -^26 -^38 ' ■ 



39 /' 



p (t3fl) _ , m (3j n (3j n (3j n (3K 

^20 - 4xix 3 m t (L; 1 3 4-^23 +-^26 +-^39) 



44x 1 x 4 m^( J Dg ) 4 £>$) 4 4x 2 x 4 m t (-Z)g ) 4 L>g } - £>g> 4 ^gio), 
where L>f } , £>^ 3) , D\f k = A , Aj , Ajfe [~f>i , &i , &2 , m t ~ ,rrig,mg, nig] . 

The F k AR ^ (k = 1 ~ 20) are written explicitly as: 

= i((x 1 x 3 m t (C 11 -C 12 ) (Ac.6) 
+x 2 x 4 m t Ci 2 - X2X 3 mpC )[-pi,pi + P2,m g ,m ii ,m il ]) 

F^ UR) = 0,(i = 3,4,...20). 
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Figure Captions 

Fig.l Feynman diagrams at the tree-level and one-loop level in the SUSY QCD for the 
gg — * ti subprocess. Fig.l (a): Tree level diagrams. Fig.l (b.l): Self-energy diagrams (for 
top-quark and gluon). Fig.l (b.2): Vertex diagrams (including tri-gluon and gluon-top-top 
interactions). Fig.l (b.3): Box diagrams (only t-channel). Dashed lines represent fi,ta in 
Fig.l (b). 

Fig. 2 (a) relative corrections to polarized and unpolarized cross sections of the ti produc- 
tion process in pp colliders as a function of y/s with input structure functions of Brodsky 
et al. U(LO); 

(b) relative corrections to polarized and unpolarized cross sections of the ti production 
process in pp colliders as a function of ^/s with input structure functions of Gliick et al. 
|]] [1C] p4| (NLO), in both above figures, solid line for the MSSM QCD correction with un- 



polarized protons, dashed line for the MSSM QCD correction with proton{+)proton{+) 
polarization, dotted line for the MSSM QCD correction with proton(—)proton(—) po- 
larization and dot-dashed line for the MSSM QCD correction with proton{+)proton{—) 
polarization; 

(c) the CP-violating parameter £cp as a function of \fs, solid line for input structure 
functions of Gliick et al(NLO), dashed line for input structure functions of Brodsky et 
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al.(LO) 

m- g = 200 GeV, m h = 250 GeV, m h = 450 GeV and 6 = <j> = 45°. 

Fig.3 (a) relative corrections to the cross section of the ti production subprocess, 5±± as a 
function of V~i, solid line for the MSSM QCD correction with gluon(+)gluon(+) polariza- 
tion and dashed line for the MSSM QCD correction with gluon{—)gluon{—) polarization, 
(b) the CP-violating parameter £cp of the subprocess as a function of VI. (c) relative 
corrections to the cross section of the subprocess J+_ as a function of y/i. 
m-g = 200 GeV, m~ tl = 250 GeV, m h = 450 GeV and 9 = <f> = 45°. 

Fig.4 (a) cross section of the ti production subprocess via gg fusion, a±± as a function 
of m-g, solid line for the MSSM QCD correction with gluon{+)gluon{+) polarization and 
dashed line for the MSSM QCD correction with gluon{—)gluon{—) polarization, (b) the 
CP-violating parameter £cp of the subprocess as a function of mg. 
m h = 100 GeV , m t ~ 2 = 450 GeV, yfs = 500 GeV and 9 = <f> = 45°. 

Fig.5 (a) relative corrections to the cross section of the ti production subprocess via 
gg fusion, d±± as a function of m^, solid line for the MSSM QCD correction with 
gluon{+)gluon(+) polarization and dashed line for the MSSM QCD correction with gluon{—)gluon{—) 
polarization, (b) the CP-violating parameter t^cp of the subprocess as a function of m^. 
m-g = 200 GeV, m h = 450 GeV, v 7 ! = 500 GeV and 9 = <p = 45°. 

Fig.6 (a) relative corrections to the cross section of the ti production subprocess via 
gg fusion, S±± as a function of m^ 2 , solid line for the MSSM QCD correction with 
gluon{+)gluon(+) polarization and dashed line for the MSSM QCD correction with gluon{—)gluon{—) 
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polarization, (b) the CP-violating parameter £cp of the subprocess as a function of mf . 
m-g = 200 GeV, m~ tl = 100 GeV and v 7 ! = 500 GeV and 6 = (j) = 45°. 
Fig. 7 (a) relative corrections to the cross section of the ti production subprocess via gg fu- 
sion, <5±± as a function of 0, solid line for the MSSM QCD correction with gluon{+)gluon(+) 
polarization and dashed line for the MSSM QCD correction with gluon{—)gluon(—) po- 
larization, (b) the CP-violating parameter ^cp of the subprocess as a function of cf). 
m-g = 200GeV,m il = 150 GeV, m h = 450 GeV ,\/| = 500 GeV and 6 = 45°. 



31 



e$(k 2 ) 
n 



t(pl) 
t(p2) 




Fig.l (a) 



/ \ 
t ' N t 




nm V(np 



n 



Fig.l (b-l) 




Fig.l (b-2) 



Fig.2(a) 



0.042 



0.037 



0.032 



0.027 



0.022 



0.017 



0.012 




1000 3000 5000 7000 9000 11000 13000 

Js (GeV) 



Fig.2(b) 



0.045 



0.04 



0.035 



0.03 



0.025 



0.02 



0.015 





unpol 




\ 

\ 

— \ \ 
\ \ 


+ + 




\ \ 

\ 

\ \ 
\ \ 
- \ \ 
\ \ 


+ - 




— \ . \ 

\ * s " 

- \ \ N 

\ \ '■• v X 
\ \ 

\ \ 
\ \. 

\ 

\ ~~- 

X 

I I I I I I I I I I I I I I I 


i i i i i i i i i i 



2000 4000 6000 8000 100001200014000 

Js (GeV) 



Fig.2(c) 



o 



0.0011 
0.001 
0.0009 
0.0008 
0.0007 
0.0006 
0.0005 
0.0004 
0.0003 
0.0002 




2000 4000 6000 8000 100001200014000 

Js (GeV) 




Fig.l (b-3) 



Fig.3(a) 





U.U A 




0.06 




0.05 




0.04 


CO 




n 




-t— ■ 


0.03 


A 




D) 


0.02 


U) 






0.01 









-0.01 




-0.02 




350 550 750 950 1150 1350 1550 1750 

Js (GeV) 



0.006 



Fig.3(b) 



0.005 



0.004 



CO 
_Q 



- 0.003 



A 



§) 0.002 



□_ 

o 



0.001 







-0.001 




350 550 750 950 1150 1350 1550 1750 

Js (GeV) 



0.03 



Fig.3(c) 



CO 
_Q 



A 
CD 

to 



0.025 



0.02 



0.015 



0.01 



0.005 



J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

350 550 750 950 1150 1350 1550 1750 

Js (GeV) 



Fig.4(a) 




100 150 200 250 300 350 400 450 500 
mass of gluino (GeV) 



0.035 
0.03 



Fig.5(a) 



0.025 



0.02 



0.015 



0.01 



0.005 




100 150 200 250 300 350 400 450 
mass of stopl (GeV) 



Fig.5(b) 

0.014 r 




100 150 200 250 300 350 400 450 



mass of stopl (GeV) 



0.035 
0.03 



Fig.6(a) 



0.025 



0.02 



0.015 



0.01 



0.005 




100 150 200 250 300 350 400 
mass of stop2 (GeV) 



Fig.6(b) 

0.014 r 




100 150 200 250 300 350 400 450 



mass of stop2 (GeV) 



0.04 



Fig.7(a) 



0.035 



0.03 



CO 

£ 0.025 

H — 1 

A 

6) 0.02 

CD 



0.015 



0.01 



0.005 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



Fig.7(b) 

0.01 r- 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 





